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1. INTRODUCTION 



The standard model [JlJ— |3J] of big bang nucle- 
osynthesis (BBN) is based on the inclusion of 
an extended nuclear network into a homogeneous 
and isotropic FRLW cosmology. There is now 
sufficient data to define the standard model in 
terms of a single parameter, namely the baryon- 
to-photon ratio, r\. Other factors, such as the 
uncertainties in reaction rates, and the neutron 
mean-life can be treated by standard statistical 
and Monte Carlo techniques to make predic- 
tions (with specified uncertainties) of the abun- 
dances of the light elements, D, 3 He, 4 He, and 7 Li. 
Even the number of neutrino flavors, N v , which 
has long been treated as a parameter can simply 
be set (=3) in defining the standard model. 

In this review, I will compare the predictions of 
BBN with the available observational determina- 
tions of the light element abundances and test for 
concordance. I will also discuss the implications 
of these results on the Galactic chemical evolu- 
tion of the light elements and on limits to particle 
properties. 

1.1. Historical Aside 

It is important to bear in mind that there has 
always been an intimate connection between BBN 
and the microwave background as a key test to 
the standard big bang model. Indeed, it was the 
formulation of BBN which predicted the existence 



of the microwave background radiation |7|. The 
argument is rather simple. BBN requires tem- 
peratures greater than 100 keV, which according 
to the standard model time-temperature relation, 
ts^MeV = 2.4/ vW, where N is the number of 
relativistic degrees of freedom at temperature T, 
corresponds to timescales less than about 200 s. 
The typical cross section for the first link in the 
nuclcosynthetic chain is 
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D + 7) ~ 5 x 10 _20 cm 3 /s 
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This implies that it was necessary to achieve a 
density 



n 10 17 citT 3 
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(2) 



The density in baryons today is known approx- 
imately from the density of visible matter to be 
riBo ~ 10 -7 cm~ 3 and since we know that that 
the density n scales as R~ 3 ~ T 3 , the tempera- 
ture today must be 



T = (n Bo /n) 1/3 T } 
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thus linking two of the most important tests of 
the Big Bang theory. 

2. THEORY 

Conditions for the synthesis of the light ele- 
ments were attained in the early Universe at tem- 
peratures T < 1 MeV. At somewhat higher tern- 
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peratures, weak interaction rates were in equilib- 
rium. In particular, the processes 
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fix the ratio of number densities of neutrons to 
protons. At T > 1 MeV, (n/p) ~ 1. 

As the temperature fell and approached the 
point where the weak interaction rates were no 
longer fast enough to maintain equilibrium, the 
neutron to proton ratio was given approximately 
by the Boltzmann factor, (n/p) ~ e~ Am / T , where 
Am is the neutron-proton mass difference. The 
final abundance of He is very sensitive to the 
(n/p) ratio 
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Freeze out occurs at slightly less than an MeV 
resulting in (n/p) ~ 1/6 at this time. 

The nucleosynthesis chain begins with the for- 
mation of Deuterium through the process, p+n — > 
D +7. However, because the large number of pho- 
tons relative to nuclcons, r\~ x = n 1 /uB ~ 10 10 , 
Deuterium production is delayed past the point 
where the temperature has fallen below the Deu- 
terium binding energy, Eb = 2.2 MeV (the aver- 
age photon energy in a blackbody is E 1 ~ 2. IT). 
This is because there are many photons in the 
exponential tail of the photon energy distribu- 
tion with energies E > Eb despite the fact that 
the temperature or E 1 is less than Eb- During 
this delay, the neutron-to-proton ratio drops to 
(n/p) ~ 1/7. 

The dominant product of big bang nucleosyn- 
thesis is 4 Hc resulting in an abundance of close to 
25% by mass. Lesser amounts of the other light 
elements are produced: D and 3 He at the level 
of about 10~ 5 by number, and 7 Li at the level of 
10 -10 by number. The resulting abundances of 
the light elements are shown in Figure [3], which 
concentrate on the range in 7710 between 1 and 10. 
The curves for the 4 He mass fraction, Y, bracket 
the computed range based on the uncertainty of 
the neutron mean-life which has been taken as 
t„ = 887 ± 2 s. Uncertainties in the produced 7 Li 



abundances have been adopted from the results in 
Hata et al. [^). Uncertainties in D and 3 He pro- 
duction are small on the scale of this figure. The 
boxes correspond to the observed abundances and 
will be discussed below. 



D, 3He 



'Li 




Figure 1. The light element abundances from big 
bang nucleosynthesis as a function of 7710 = 10 10 ?7. 



3. Data 
3.1. 4 He 

The primordial 4 He abundance is best deter- 
mined from observations of Hell — > Hel recombi- 
nation lines in extragalactic HII (ionized hydro- 
gen) regions. There is now a good collection of 
abundance information on the 4 He mass fraction, 
Y, O/H, and N/H in over 70 |-[l0) such regions. 
Since He is produced in stars along with heav- 
ier elements such as Oxygen, it is then expected 
that the primordial abundance of 4 He can be de- 
termined from the intercept of the correlation be- 
tween Y and O/H, namely Y p = V(0/H -> 0). A 
detailed analysis of the data including that in JTof ] 
found an intercept corresponding to a primordial 
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abundance Y p = 0.234 ± 0.002 ± 0.005 [|T). This 
was updated to include the most recent results 
of [[l2| in [^3|. The result (which is used in the 
discussion below) is 



Y p = 0.238 ± 0.002 ±0.005 



(5) 



The first uncertainty is purely statistical and the 
second uncertainty is an estimate of the system- 
atic uncertainty in the primordial abundance de- 
termination [jy]. The solid box for 4 He in Fig- 
ure [l] represents the range (at 2er sta t) from (||). 
The dashed box extends this by including the 
systematic uncertainty. A somewhat lower pri- 
mordial abundance of Y p — 0.235 ± .003 is found 
by restricting to the 36 most metal poor regions 
|l3[ . These results are consistent with those of 
a Bayesian analysis [[l4j based on the 32 points 
of lowest metallicity. These have been used to 
calculate a likelihood function for which the peak 
occurs at Y p = 0.238 and the most likely spread is 
w = 0.009. The 95% CL upper limit to Y p in this 
case is 0.245. For further details on this approach 
see R. 



24%. In particular, even in systems in which an 
element such as Oxygen, which traces stellar ac- 
tivity, is observed at extremely low values (com- 
pared with the solar value of O/H « 8.5 x 10~ 4 ), 
the 4 He abundance is nearly constant. This is far 
different from all other element abundances (with 
the exception of 7 Li as we will see below). For ex- 
ample, in Figure [| the N/H vs. O/H correlation 
is shown. As one can clearly see, the abundance 
of N/H goes to 0, as O/H goes to 0, indicating a 
stellar source for Nitrogen. 




» 8 » 4 He is Primordial! 



10* O/H 



Figure 2. The He mass fraction as determined in 
extragalactic H II regions as a function of O/H. 



A global view of the 4 He data is shown in Fig- 
ure ||. What should be absolutely apparent from 
this figure is the primordial nature of 4 He. There 
are no observations to date which indicate a 4 He 
abundance which is significantly below 23% to 



Figure 3. The Nitrogen and Oxygen abundances 
in the same extragalactic HII regions with ob- 
served 4 He shown in Figure 0. 



A more useful plot of the 4 He data is shown 
in Figure [|. Here one sees the correlation of 4 He 
with O/H and the linear regression which leads 
to primordial abundance given in Eq. (^). 

3.2. 7 Li 

The abundance of 7 Li has been determined by 
observations of over 100 hot, population-II stars, 
and is found to have a very nearly uniform abun- 
dance |l5||. For stars with a surface temperature 
T > 5500 K and a metallicity less than about 
l/20th solar (so that effects such as stellar con- 
vection may not be important), the abundances 
show little or no dispersion beyond that which is 
consistent with the errors of individual measure- 
ments. Indeed, as detailed in ref. |L6 17j, much 
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Figure 4. The Helium (Y) and Oxygen (O/H) 
abundances in extragalactic HII regions, from 
refs. an d from ref. fj"^| . Lines connect 

the same regions observed by different groups. 
The regression shown leads to the primordial 4 He 
abundance given in Eq. (||) . 




Figure 5. The Li abundance in halo stars with 
[Fc/H] < -f.3, as a function of surface temper- 
ature. The dashed line shows the value of the 
weighted mean of the plateau data. 



of the work concerning 7 Li has to do with the 
presence or absence of dispersion and whether or 
not there is in fact some tiny slope to a [Li] = log 
7 Li/H + 12 vs. T or [Li] vs. [Fe/H] relationship 
([Fe/H] is the log of the Fe/H ratio relative to the 
solar value). 

When the Li data from stars with [Fe/H] < 
-1.3 is plotted as a function of surface tempera- 
ture, one sees a plateau emerging for T > 5500 
K as shown in Figure || for the data taken from 
ref. pq| . As one can see from the figure, at high 
temperatures, where the convection zone does not 
go deep below the surface, the Li abundance is 
uniform. At lower surface temperatures, the sur- 
face abundance of Li is depleted as Li is dragged 
through the hotter interior of the star and de- 
stroyed. The lack of dispersion in the plateau 
region is evidence that this abundance is indeed 
primordial (or at least very close to it). 

I will use the value given in ref. |l(| as the 
best estimate for the mean 7 Li abundance and its 
statistical uncertainty in halo stars 

Li/H = (1.6 ± 0.1) x 10~ 10 (6) 

The small error is is statistical and is due to the 
large number of stars in which 7 Li has been ob- 
served. The solid box for 7 Li in Figure II] repre- 



sents the 2cr st at range from . There is however 
an important source of systematic error due to 
the possibility that Li has been depleted in these 
stars, though the lack of dispersion in the Li data 
limits the amount of depletion. In addition, stan- 
dard stellar models jll| predict that any depletion 
of Li would be accompanied by a very severe de- 
pletion of 6 Li. Until recently, 6 Li had never been 
observed in hot pop II stars. The observation ]Tg|] 
of 6 Li (which turns out to be consistent with its 
origin in cosmic-ray nucleosynthesis) is another 
good indication that 7 Li has not been destroyed 
in these stars pp|-|22]] . 

Aside from the big bang, Li is produced to- 
gether with Be and B in accelerated particle in- 
teractions such as cosmic ray spallation of C,N,0 
by protons and a-particles. Li is also produced 
by a — a fusion. Be and B have been observed 
in these same pop II stars and in particular there 
are a dozen or so stars in which both Be and 
7 Li have been observed. Thus Be (and B though 
there is still a paucity of data) can be used as a 
consistency check on primordial Li Based 
on the Be abundance found in these stars, one 
can conclude that no more than 10-20% of the 
7 Li is due to cosmic ray nucleosynthesis leaving 
the remainder (an abundance near 10 -10 ) as pri- 
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mordial. A similar conclusion was reached in p7| . 
The dashed box in Figure [l], accounts for the pos- 
sibility that as much as half of the primordial Li 
has been destroyed in stars, and that as much 
as 20% of the observed 7 Li may have been pro- 
duced in cosmic ray collisions rather than in the 
Big Bang. For 7 Li, the uncertainties are clearly 
dominated by systematic effects. 

4. Likelihood Analyses 

At this point, having established the primor- 
dial abundance of at least two of the light ele- 
ments, 4 He and 7 Li, with reasonable certainty, it 
is possible to test the concordance of BBN the- 
ory with observations. Monte Carlo techniques 
have proven to be a useful form of analysis in 
this regard Two elements are sufficient for 

not only constraining the one parameter theory 
of BBN, but also for testing for consistency |E3 . 
The procedure begins by establishing likelihood 
functions for the theory and observations. For 
example, for 4 He, the theoretical likelihood func- 
tion takes the form 

L BBN (Y, Ybbn) = e-^-W.) W (7) 

where 1bbn(??) is the central value for the 4 He 
mass fraction produced in the big bang as pre- 
dicted by the theory at a given value of 77. <i\ 
is the uncertainty in that value derived from the 
Monte Carlo calculations || and is a measure of 
the theoretical uncertainty in the big bang calcu- 
lation. Similarly one can write down an expres- 
sion for the observational likelihood function. As- 
suming Gaussian errors, the likelihood function 
for the observations would take a form similar to 
that in (g). 

A total likelihood function for each value of rj 
is derived by convolving the theoretical and ob- 
servational distributions, which for 4 He is given 
by 

= J dYL BBN (Y, y BB N (v)) Lo(Y, Y ) (8) 

An analogous calculation is performed |2j] for 
7 Li. The resulting likelihood functions from the 
observed abundances given in Eqs. (|) and (§) 
is shown in Figure §. As one can see there is 



very good agreement between 4 He and 7 Li in the 
range of 7710 — 1-5 - 5.0. The double peaked na- 
ture of the 7 Li likelihood function is due to the 
presence of a minimum in the predicted lithium 
abundance. For a given observed value of 7 Li, 
there are two likely values of 77. 
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Figure 6. Likelihood distribution for each of 4 He 
and 7 Li, shown as a function of 77. The one- 
peak structure of the 4 He curve corresponds to 
the monotonic increase of Y p with 77, while the 
two peaks for 7 Li arise from the minimum in the 
7 Li abundance prediction. 




Figure 7. Combined likelihood for simultaneously 
fitting 4 He and 7 Li, as a function of 77. 



The combined likelihood, for fitting both ele- 
ments simultaneously, is given by the product of 



G 



the two functions in Figure |6j and is shown in 
Figure 0. The 95% CL region covers the range 
1.55 < 7710 < 4.45, with the two peaks occurring 
at 7710 =1.9 and 3.5. This range corresponds to 
values of Qb between 

0.006 < n B h 2 < .016 (9) 

5. More Data 

5.1. D and 3 He 

Because there are no known astrophysical sites 
for the production of Deuterium, all observed D 
is assumed to be primordial. As a result, any firm 
determination of a Deuterium abundance estab- 
lishes an upper bound on r\ which is robust. 

Deuterium abundance information is available 
from several astrophysical environments, each 
corresponding to a different evolutionary epoch. 
In the ISM, corresponding to the present epoch, 
an often quoted measurement for D/H is |2f| 

(D/H)i SM = 1.60 ± 0.09±Ho x !0~ 5 (10) 

This measurement allow us to set the upper limit 
to 7710 < 9 and is shown by the lower right of the 
solid box in Figure 0. There are however, serious 
questions regarding homogeneity of this value in 
the ISM. There may be evidence for considerable 
dispersion in D/H |p6| , [27f as is the case with 3 He 
& 

The solar abundance of D/H is inferred from 
two distinct measurements of 3 He. The solar 
wind measurements of 3 He as well as the low tem- 
perature components of step-wise heating mea- 
surements of 3 He in meteorites yield the preso- 
lar (D + 3 He)/H ratio, as D was efficiently 
burned to 3 He in the Sun's pre-main-sequence 
phase. These measurements indicate that |29, 
|FJ ((D+ 3 He)/H) Q = (4.1 ±0.6 ±1.4) x lCT^. 
The high temperature components in meteorites 
are believed to yield the true solar 3 He/H ratio 
of HQ ( 3 He/H) Q = (1.5 ± 0.2 ± 0.3) x 10" 5 . 
The difference between these two abundances re- 
veals the presolar D/H ratio, giving, 

(D/H)q w (2.6 ± 0.6 ± 1.4) x 10~ 5 (11) 

This value for presolar D/H is consistent with 
measurements of surface abundances of HD on 
Jupiter D/H = 2.7 ± 0.7 x lO" 5 ||. 



Finally, there have been several reported mea- 
surements of D /H in high redshift quasar absorp- 
tion systems. Such measurements are in princi- 
ple capable of determining the primordial value 
for D/H and hence 77, because of the strong and 
monotonic dependence of D/H on 77. However, at 
present, detections of D/H using quasar absorp- 
tion systems do not yield a conclusive value for 
D/H. As such, it should be cautioned that these 
values may not turn out to represent the true pri- 
mordial value and it is very unlikely that both arc 
primordial and indicate an inhomogeneity |32[] (a 
large scale inhomogeneity of the magnitude re- 
quired to placate all observations is excluded by 
the isotropy of the microwave background radi- 
ation). The first of these measurements |33| in- 
dicated a rather high D/H ratio, D/H w 1.9 - 
2.5 xlO -4 . Other high D/H ratios were reported 
in |34|. More recently, a similarly high value of 
D/H = 2.0 ±0.5 x 10~ 4 was reported in a rel- 
atively low redshift system (making it less sus- 
pect to interloper problems) |55|. This was con- 
firmed in |36| where a 95% CL lower bound to 
D/H was reported as 8 xlO -5 . However, there 
are reported low values of D/H in other such sys- 
tems |37j with values of D/H originally reported 
as low as ~ 2.5 x 10~ 5 , significantly lower than 
the ones quoted above. The abundance in these 
systems has been revised upwards to about 3.4 
±0.3 x 10" 5 |||. However, it was also noted ||| 
that when using mesoturbulent models to account 
for the velocity field structure in these systems, 
the abundance may be somewhat higher (3.5 - 
5 xl0~ 5 ). This may be quite significant, since 
at the upper end of this range (5 xlO -5 ) all of 
the element abundances are consistent as will be 
discussed shortly. I will not enter into the de- 
bate as to which if any of these observations may 
be a better representation of the true primordial 
D/H ratio. The status of these observations are 
more fully reviewed in |27| . The upper range of 
quasar absorber D /H is shown by the dashed box 
in Figure ^. 

There are also several types of 3 He measure- 
ments. As noted above, meteoritic extractions 
yield a presolar value for 3 He/H. In addition, 
there are several ISM measurements of 3 He in 
galactic HII regions |28) which show a wide dis- 
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persion which may be indicative of pollution or 
a bias @ ( 3 He/H) fflI ~ 1 - 5 x 10~ 5 . There 
is also a recent ISM measurement of 3 He jll) 
with ( 3 He/H) ISM = 2.11 j x 10~ 5 . Finally there 
are observations of 3 He in planetary nebulae Jl2) 
which show a very high 3 He abundance of 3 Hc/H 
~ 10~ 3 . None of these determinations represent 
the primordial 3 He abundance, and as will be 
discussed below, their relation to the primordial 
abundance is heavily dependent on both stellar 
and chemical evolution. 

6. More Analysis 

It is interesting to compare the results from 
the likelihood functions of 4 He and 7 Li with that 
of D/H. Since D and 3 He are monotonic func- 
tions of rj, a prediction for rj, based on 4 He and 
7 Li, can be turned into a prediction for D and 
3 He. The corresponding 95% CL ranges are D/H 
= (4.1-25) xlCT 5 and 3 He/H= (1.2-2.6) XHT 5 . 
If we did have full confidence in the measured 
value of D/H in quasar absorption systems, then 
we could perform the same statistical analysis us- 
ing 4 He, 7 Li, and D. To include D/H, one would 
proceed in much the same way as with the other 
two light elements. We compute likelihood func- 
tions for the BBN predictions as in Eq. (0) 
and the likelihood function for the observations. 
These are then convolved as in Eq. (^). 

Using D/H = (2.0 ±0.5) x 10~ 4 as indicated in 
the high D/H systems, we can plot the three like- 
lihood functions including L^ )tal (^) in Figure |[ 
It is indeed startling how the three peaks, for D, 
4 He and 7 Li are in excellent agreement with each 
other. In Figure ^, the combined distribution is 
shown. We now have a very clean distribution 
and prediction for 77, rj W — 1.8^ 3 ' correspond- 
ing to tt B h 2 = .007l;°2i- Tne absence of any 
overlap with the high-?7 peak of the 7 Li distribu- 
tion has considerably lowered the upper limit to 
77. Overall, the concordance limits in this case are 
dominated by the Deuterium likelihood function. 

If instead, we assume that the low value [p8| 
of D/H = (3.4 ± 0.3) x 10" 5 is the primordial 
abundance, then we can again compare the like- 
lihood distributions as in Figure ||, now substi- 
tuting the low D/H value. As one can see from 
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Figure 8. As in Figure ^, with the addition of the 
likelihood distribution for D/H assuming "high" 
D/H. 
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Figure 9. Combined likelihood for simultaneously 
fitting 4 He and 7 Li, and D as a function of ry from 
Figure ^. 
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Figure [l^, there is now hardly any overlap be- 
tween the D and the 7 Li and 4 He distributions. 
The combined distribution shown in Figure [ll] is 
compared with that in Figure ^|. In this case, 
D/H is just compatible (at the 2 a level) with the 
other light elements, and the peak of the likeli- 
hood function occurs at 7710 = 4.8i ; 6 . Though 
one can not use this likelihood analysis to prove 
the correctness of the high D /H measurements or 
the incorrectness of the low D/H measurements, 
the analysis clearly shows the difference in com- 
patibility between the two values of D /H and the 
observational determinations of 4 He and 7 Li. To 
make the low D /H measurement compatible, one 
would have to argue for a shift upwards in 4 He 
to a primordial value of 0.247 (a shift by 0.009) 
which is not warranted at this time by the data, 
and a 7 Li depletion factor of about 2, which is 
close to recent upper limits to the amount of de- 
pletion g§|2l]]. 
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Figure 10. As in Figure |], with the likelihood 
distribution for low D/H. 



It is important to recall however, that the true 
uncertainty in the low D/H systems might be 
somewhat larger. If we allow D/H to be as large 
as 5 x 10 -5 , the peak of the D/H likelihood func- 
tion shifts down to 7/10 ~ 4. In this case, there 
would be a near perfect overlap with the high 77 
7 Li peak and since the 4 He distribution function 
is very broad, this would be a highly compatible 
solution. 




Figure 11. Combined likelihood for simultane- 
ously fitting 4 He and 7 Li, and low D/H as a func- 
tion of 77. The dashed curve represents the com- 
bined distribution shown in Figure ^. 



7. Chemical Evolution 

Because we can not directly measure the pri- 
mordial abundances of any of the light element 
isotopes, we are required to make some assump- 
tions concerning the evolution of these isotopes. 
As has been discussed above, 4 He is produced 
in stars along with Oxygen and Nitrogen. 7 Li 
can be destroyed in stars and produced in several 
(though still uncertain) environments. D is to- 
tally destroyed in the star formation process and 
3 He is both produced and destroyed in stars with 
fairly uncertain yields. It is therefore preferable, 
if possible to observe the light element isotopes in 
a low metallicity environment. Such is the case 
with 4 He and 7 Li. These elements are observed 
in environments which are as low as l/50th and 
1/ 1000th solar metallicity respectively and we can 
be fairly assured that the abundance determina- 
tions of these isotopes are close to primordial. If 
the quasar absorption system measurements of 
D/H stabilize, then this too may be very close to 
a primordial measurement. Otherwise, to match 
the solar and present abundances of D and 3 He to 
their primordial values requires a model of galac- 
tic chemical evolution. 

The main inputs to chemical evolution models 
are: 1) The initial mass function, 4>(m), indicat- 
ing the distribution of stellar masses. Typically, 
a simple power law form for the IMF is chosen, 
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4>{m) ~ m~ x , with x ~ —2.7. This is a fairly good 
representation of the observed distribution, par- 
ticularly at larger masses. 2) The star formation 
rate, -0- Typical choices for a SFR are ip(t) oc cr or 
<7 2 or even a straight exponential e~ t l r . a is the 
fraction of mass in gas, M gas /M to t- 3) The pres- 
ence of infalling or outflowing gas; and of course 
4) the stellar yields. It is the latter, particularly 
in the case of 3 He, that is the cause for so much 
uncertainty. Chemical evolution models simply 
set up a series of evolution equations which trace 
desired quantities. 

Deuterium is always a monotonically decreas- 
ing function of time in chemical evolution models. 
The degree to which D is destroyed, is however 
a model dependent question which depends sen- 
sitively on the IMF and SFR. The evolution of 
3 He is however considerably more complicated. 
Stellar models predict that substantial amounts 
of 3 He are produced in stars between 1 and 3 
Mq. For example, in the models of Iben and Tru- 
ran j44j a 1 M Q star will yield a 3 He abundance 
which is nearly three times as large as its initial 
(D+ 3 He) abundance. It should be emphasized 
that this prediction is in fact consistent with the 
observation of high 3 He/H in planetary nebulae 

E§. 

However, the implementation of standard 
model 3 He yields in chemical evolution models 
leads to an overproduction of 3 He/H particularly 
at the solar epoch jlO 45 1 . While the overproduc- 
tion is problematic for any initial value of D/H, 
it is particularly bad in models with a high pri- 
mordial D/H. In Scully et al. [Q, a dynamically 
generated supernovae wind model was coupled to 
models of galactic chemical evolution with the 
aim of reducing a primordial D/H abundance of 
2 xlO -4 to the present ISM value without over- 
producing heavy elements and remaining consis- 
tent with the other observational constraints typ- 
ically imposed on such models. In Figure [l^, the 
evolution of D/H and 3 He/H is shown as a func- 
tion of time in one of the models with signifi- 
cant Deuterium destruction factors (see ref [[Ig| 
for details). While such a model can success- 
fully account for the evolution of D /H (and other 
standard chemical evolution tracers), as one can 
plainly see, 3 He is grossly overproduced (the Deu- 



terium data is represented by squares and 3 He by 
circles). 
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Figure 12. The evolution of D/H and 3 He/H 
with time in units of 10~ J . The assumed primor- 
dial abundance of D/H is 2 xlO~ 4 . The dashed 
3 Hc curve shows the evolution with standard 3 He 
yields, while the solid curve shows the effect of 
the reduced yields of [48[. 



The overproduction of 3 He relative to the so- 
lar meteoritic value seems to be a generic feature 
of chemical evolution models when 3 He produc- 
tion in low mass stars is included. This result ap- 
pears to be independent of the chemical evolution 
model and is directly related to the assumed stel- 
lar yields of 3 He. It has recently been suggested 
that at least some low mass stars may indeed be 
net destroyers of 3 He if one includes the effects of 
extra mixing below the conventional convection 
zone in low mass stars on the red giant branch 
|47||4j|. The extra mixing does not take place 
for stars which do not undergo a Helium core 
flash (i.e. stars > 1.7 - 2 M Q ). Thus stars with 
masses less than 1.7 M are responsible for the 
3 He destruction. Using the yields of Boothroyd 
and Malaney |48[j , it was shown ji^] that these 
reduced 3 He yields in low mass stars can account 
for the relatively low solar and present day 3 He/H 
abundances observed. In fact, in some cases, 3 He 
was underproduced. To account for the 3 He evo- 
lution and the fact that some low mass stars must 
be producers of 3 He as indicated by the planetary 
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nebulae data, it was suggested that the new yields 
apply only to a fraction (albeit large) of low mass 
stars |l9],[50| . The corresponding evolution of 
D/H and 3 He/H using the reduced yields is shown 
in Figure |lj. 

The models of chemical evolution discussed 
above indicate that it is possible to destroy signifi- 
cant amounts of Deuterium and remain consistent 
with chemical evolutionary constraints. To do so 
however, comes with a price. Large Deuterium 
destruction factors require substantial amounts 
of stellar processing, which at the same time pro- 
duce heavy elements. To keep the heavy element 
abundances in the Galaxy in check, significant 
Galactic winds enriched in heavy elements must 
be incorporated. In fact there is some evidence 
that enriched winds were operative in the early 
Galaxy. In the X-ray cluster satellites observed 
by Mushotzky et al. jsij and Loewenstein and 
Mushotzky |32j the mean Oxygen abundance was 
found to be roughly half solar. This corresponds 
to a near solar abundance of heavy elements in 
the inter-Galactic medium, where apparently lit- 
tle or no star formation has taken place. 

If our Galaxy is typical in the Universe, then 
the models of the type discussed above would in- 
dicate that the luminosity density of the Universe 
at high redshift should also be substantial aug- 
mented relative to the present. Recent observa- 
tions of the luminosity density at high redshift 
]53fl are making it possible for the first time to 
test models of cosmic chemical evolution. The 
high redshift observations, are very discrimina- 
tory with respect to a given SFR |54[] . Models in 
which the star formation rate is proportional to 
the gas mass fraction (these are common place in 
Galactic chemical evolution) have difficulties to 
fit the multi-color data from z = to 1. This in- 
cludes many of the successful Galactic infall mod- 
els. In contrast, models with a steeply decreas- 
ing SFR are favored. In Figure [l3|, the predicted 
luminosity density based on the model with evo- 
lution shown in Figure ^ from |4^| , as compared 
with the observations (see ref. |5J] for details). 

While it would be premature to conclude that 
all models with large Deuterium destruction fac- 
tors are favored, it does seem that models which 
do fit the high redshift data destroy significant 




Figure 13. The tricolor luminosity densities (UV, 
B and IR) at A = 0.28,0.44 and 1.0 fim, in units 
of (h/.5) WHz -1 Mpc -3 as a function of redshift 
for a model shown in |l2| which destroys significant 
amounts of D/H. The data are taken from f53| . 
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amounts of D/H. On the other hand, we can 
not exclude models which destroy only a small 
amount of D/H as Galactic models of chemical 
evolution. In this case, however the evolution 
of our Galaxy is anomalous with respect to the 
cosmic average. If the low D/H measurements 
|37|]38j 1 hold up, then it would seem that our 
Galaxy also has an anomalously high D/H abun- 
dance. That is we would predict in this case that 
the present cosmic abundance of D/H is signifi- 
cantly lower than the observed ISM value. (This 
conclusion assumes that the ISM D /H abundance 
is representative of the present epoch. If new re- 
sults [^6| which show values as low as a few x 10 ~ 6 
are more typical, then low primordial D/H would 
fit the high redshift data and models with a high 
degree of D destruction much better.) If the high 
D/H observations p3[-p5[ hold up, we could con- 
clude that our Galaxy is indeed representative of 
the cosmic star formation history. 

8. Constraints from BBN 

Limits on particle physics beyond the standard 
model are mostly sensitive to the bounds imposed 
on the 4 He abundance. As is well known, the 4 He 
abundance is predominantly determined by the 
neutron-to-proton ratio just prior to nucleosyn- 
thesis and is easily estimated assuming that all 
neutrons are incorporated into 4 He. As discussed 
earlier, the neutron-to-proton ratio is fixed by its 
equilibrium value at the freeze-out of the weak 
interaction rates at a temperature Tf ~ 1 MeV 
modulo the occasional free neutron decay. Fur- 
thermore, freeze-out is determined by the compe- 
tition between the weak interaction rates and the 
expansion rate of the Universe 



G F 2 T f 5 ~ r wk (2» = H(Tf) ~ y/G N NT f 2 (12) 

where N counts the total (equivalent) number of 
relativistic particle species. The presence of ad- 
ditional neutrino flavors (or any other relativistic 
species) at the time of nucleosynthesis increases 
the overall energy density of the Universe and 
hence the expansion rate leading to a larger value 
of Tf, (n/p), and ultimately Y p . Because of the 
form of Eq. (112) it is clear that just as one can 
place limits [pa] on N, any changes in the weak or 



gravitational coupling constants can be similarly 
constrained (for a recent discussion see ref. |pql ). 

As discussed above, the limit on N v comes 
about via the change in the expansion rate given 
by the Hubble parameter, 



H 2 



[N SM + t^A^T 4 



-P = 



3 90 L "'° 1V1 ' 8~""'' (13) 

when compared to the weak interaction rates. 
Here Asm refers to the standard model value for 
N. At T ~ 1 MeV, A SM = 43/4. Additional 
degrees of freedom will lead to an increase in 
the freeze-out temperature eventually leading to 
a higher 4 He abundance. In fact, one can param- 
eterize the dependence of Y on N v by 

Y = 0.2262 + 0.0131(A„ - 3) + 0.0135 In 7710 (14) 

in the vicinity of 7710 ~ 2. Eq. ( |l4| ) also shows 
the weak (log) dependence on 77. However, rather 
than use fl4| ) to obtain a limit, it is preferable to 
use the likelihood method. 

Just as 4 He and 7 Li were sufficient to determine 
a value for 77, a limit on N v can be obtained as 
well [|4K5|j5|,§. The likelihood approach utilized 
above can be extended to include N v as a free 
parameter. Since the light clement abundances 
can be computed as functions of both 77 and N v , 
the likelihood function can be defined by J 57 re- 
placing the quantity Ybbn (v) m e( l- 

©with 

Ybbn(v,N u ) to obtain L Hc to tai(?7, N„). Again, 
similar expressions are needed for 7 Li and D. 

A three-dimensional view of the combined like- 
lihood functions |58j is shown in Figure [TJ. In 
this case the high and low 77 maxima of Figure [7], 
show up as peaks in the L — rj — N u space. The 
likelihood function is labeled L47 (and L247 when 
D/H is included). In Figures [FJ and [L6| the corre- 
sponding likelihood functions £247 with high and 
low D/H are shown. Once again one sees an ef- 
fect of including D /H is to eliminate one of the 7 Li 
peaks. Furthermore, unlike the case discussed in 
section 6, the likelihood distribution for low D/H 
is just as strong as that for high D/H, albeit at a 
lower value of N v . 

The peaks of the distribution as well as the 
allowed ranges of 77 and N v are more easily dis- 
cerned in the contour plots of Figures and|l| 
which show the 50%, 68% and 95% confidence 
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level contours in L47 and L247 for high and low 
D/H as indicated. The crosses show the loca- 
tion of the peaks of the likelihood functions. L47 
peaks at N u = 3.2, ij w = 1.85 and at N v = 2.6, 
7710 = 3.6. The 95% confidence level allows the 
following ranges in 77 and N v 

1.7 <N V < 4.3 1.4 < 77x0 < 4.9 (15) 

Note however that the ranges in 77 and N v are 
strongly correlated as is evident in Figure ft7j . 

With high D/H, i 2 47 peaks at = 3.3, and 
also at 7710 = 1-85. In this case the 95% contour 
gives the ranges 

2.2 <N U < AA 1.4 < 7710 < 2.4 (16) 

Note that within the 95% CL range, there is also a 
small area with 7710 = 3.2 — 3.5 and N v — 2.5 — 2.9. 

Similarly, for low D/H, L247 peaks at N u — 
2.4, and 7710 = 4.55. The 95% CL upper limit 
is now N v < 3.2, and the range for 77 is 3.9 < 
7710 < 5.4. It is important to stress that with 
the increase in the determined value of D /H |3£| 
in the low D/H systems, these abundances are 
now consistent with the standard model value of 
N v = 3 at the 2 a level. 




Figure 15. L 2 47(N u ,r]) for observed abundances 
including high D/H. 




Figure 14. L^j{N u ^ri) for observed abundances 
given by eqs. (|| and ||). 



Figure 16. L 2 47(N u ,r]) for observed abundances 
including low D/H. 



13 



5 



4 



2 



. ' ' ' 1 




D/H = (2.0+0.5)X10^z 






I 47 j 












Iftl : 


: , , , 1 







1 10 



Figure 17. 50%, 68% & 95% C.L. contours of L 47 
and L247 where observed abundances are given by 
eqs. ^ and and high D/H. 



9. Summary 

To summarize, I would assert that one can con- 
clude that the present data on the abundances of 
the light element isotopes are consistent with the 
standard model of big bang nucleosynthesis. Us- 
ing the isotopes with the best data, 4 He and 7 Li, 
it is possible to constrain the theory and obtain a 
best set of values for the baryon-to-photon ratio 
of 7710 and the corresponding value for fish 2 

1.55 < r)w < 4.45 95%CL 

.006 < n B h 2 < .016 95%CL 1 ' 

For 0.4 < h < 1, we have a range .006 < B < 
.10. This is a rather low value for the baryon den- 
sity and would suggest that much of the galactic 
dark matter is non-baryonic |30|. These predic- 
tions arc in addition consistent with recent ob- 
servations of D/H in quasar absorption systems 
which show a high value. Difficulty remains how- 
ever, in matching the solar 3 He abundance, sug- 
gesting a problem with our current understanding 
of galactic chemical evolution or the stellar evo- 
lution of low mass stars as they pertain to 3 He. 



D/H = (3.4±0.3)X10~" 

' 47 

L 247 



Figure 18. 50%, 68% & 95% C.L. contours of L47 
and L247 where observed abundances are given by 
eqs. (H and ||), and low D/H. 
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